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Abstract 

We report angle-resolved Raman scattering observations of the temperature 
dependent Landau damping of the acoustic plasmon in an electron bilayer 
system realised in a GaAs double quantum well structure. Corresponding 
calculations of the charge-density excitation spectrum of the electron bilayer 
using forms of the random phase approximation (RPA), and the static local 
field formalism of Singwi, Tosi, Land and Sjolander (STLS) extended to in- 
corporate non-zero electron temperature T c and phenomenological damping, 
are also presented. The STLS calculations include details of the temperature 
dependence of the intra- and inter-layer local field factors and pair-correlation 
functions. Good agreement between experiment and the various theories is ob- 
tained for the acoustic plasmon energy and damping for T c < If/2, the Fermi 
temperature. However, contrary to current expectations, all of the calcula- 
tions show significant departures from our experimental data for T e > Tp/2. 
From this, we go on to demonstrate unambiguously that real local field fac- 
tors fail to provide a physically accurate description of exchange correlation 
behaviour in low dimensional electron gases. Our results suggest instead that 
one must resort to a dynamical local field theory, characterised by a complex 
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field factor to provide a more accurate description. 
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I. INTRODUCTION 



Low- dimensional electron gas systems confined in semiconductor heterostructures pro- 
vide ideal systems for the study of electron-electron interactions, providing a high degree of 
structural quality and control. The usual theoretical approach adopted for the description 
of the collective excitation spectra of such systems is the Random Phase Approximation 
(RPA). However, the RPA does not include the effects of exchange and correlation, which 
are expected to be much more important in two-dimensional electron gas (2DEG) systems 
than in 3D, and to increase in importance as the number density of the system decreases.0 
To account for these effects, many authors have gone beyond the RPA for a description of 
the excitation spectra of such systems, by adopting the local field theories of Hubbard! 
and Singwi, Tosi, Land and Sjolander (STLS),@ which include corrections for exchange and 
correlation in a simple, physically motivated mannerjl! the success of these approaches can 
be seen by comparing their results against Monte Carlo calculations.0 There have also been 
reports of a number of alternative approaches for including exchange-correlation effects in the 
excitation spectra of two-dimensional electron systems by adopting: (i) a non-local approach 
within a Hartree-Fock approximation for the evaluation of the irreducible polarisability (see, 
e.g., Refs. [11] and [12] for the case of intersubband excitations in quantum wells); (ii) alterna- 
tive static local field approaches such as the quasi-localized charge methodjll (iii) a dynamic 
local field approach within the quantum STLS theory.O However, the static local field STLS 
approach has remained a firm favourite in the published literature. Tests of local field the- 
ories for three-dimensional systems have been reported for, e.g., AluminiumEl but there are 
relatively few reliable experimental tests of these widely used theories in the density range 
applicable to semiconductor heterostructures. We report here such an experimental test and 
demonstrate that static local field corrections are unable to provide an accurate description 
of exchange correlation behaviour in low dimensional electron gases. Indeed we will argue 
that our results point to the need to examine dynamical local field theories, embodied by a 
complex local field factor to achieve a more accurate description of our experimental results. 



3 



Electron bilayers (two parallel 2DEGs in close proximity, produced by MBE-growth of a 
GaAs/AlGaAs modulation-doped double quantum well heterostructure) have been identified 
as useful systems for the study of electron correlations in low-dimensional electron systems 
because the inter-layer Coulomb interaction can effectively counterbalance the kinetic energy 
of the electrons, allowing many body effects to dominate. Such a system supports two plas- 
mon modes corresponding to the in-phase (optic plasmon, OP) and out-of-phase (acoustic 
plasmon, AP) intra-sub-band oscillations of the charge densities in each layer.El Liu et al. 
have shown that electron correlations can push the acoustic plasmon curve completely into 
the single-particle excitation (SPE) spectrum, resulting in Landau damping of the mode at 
much smaller wavevectors than predicted by the RPA.i The acoustic plasmon may therefore 
be used as a probe of exchange-correlation effects in 2DEG systems, and of Landau damping 
effects. 

Raman spectroscopy is a powerful tool for the study of the electronic excitation spectra 
of low- dimensional semiconductor systems.00 For example, measurements of intersubband 
excitations in quantum wells have indicated the importance of exchange-correlation in these 
systemsJil Applying angle resolved Raman techniques to electron bilayers, we have recently 
determined the dispersion of the acoustic plasmon of a double quantum well structure-SIll 
However, the in-plane wavevectors q accessible (q < 0.15 kp [the Fermi wavevector], for 
2DEG areal densities N s ~ 2x 10 11 cm~ 2 ), determined by the wavelength of the incident light, 
are too small to probe exchange-correlation corrections to the AP dispersion, as predicted in 
Ref. |5]. Attempts to extend the wavevector range by overlaying a grating to couple in higher 
wavevector transfersil did not provide significantly clear results for any conclusions to be 
drawn; the superposition of the folded AP and OP branches in the reduced Brillouin zone 
of the grating leads to ambiguities in any analysis of the form of the plasmon dispersions^ 
However, substantial overlap between the SPE continuum and the AP can be obtained by 
increasing the electron temperature, resulting in the Landau damping of the plasmon, as we 
have demonstrated recently.^ 

In this earlier work we characterised the dependence of the Landau damping of the AP 



mode in a bilayer system on electron temperature T e by comparing the AP peak widths 
measured using Raman spectroscopy with those calculated using the RPA and the RPA 
including a Hubbard (static, zero-T e ) local field correction (H-RPA) to account for exchange- 
correlation effect s.S The H-RPA agreed better with the data for low temperatures (T c < 
Tp/2, the Fermi temperature) than did the RPA, suggesting that exchange-correlation effects 
are important even for the comparatively high 2DEG densities used; but for higher T e , 
the RPA gave a better fit, suggesting that at higher temperatures (T e > Tp/2) exchange- 
correlation effects become less significant and the RPA the more appropriate approach. 

The temperature dependence of the Coulomb drag between two 2DEGs has also been 
employed recently to study the effects of electron correlations in bilayer systems.@il The 
drag mechanism was attributed to Landau damping of the AP,il0 but although inclusion 
of a K Hubbard local field correction provided a better description of the data than that 
given by the RPA,@ overall agreement between experiment and theory was not achieved. In 
both these Raman and Coulomb drag studies0@ it was suggested that using a static local 
field correction which was Independent might be important, and the present work addresses 
this point in detail by comparing experimental data with corresponding calculations. 

We present here Raman scattering measurements of the Landau damping induced width 
of the AP as a function of T e , and hence obtain a measure of the overlap of the AP with the 
SPE continuum and of the importance of exchange and correlation in these systems. The 
data are modelled using calculations (taking full account of the finite thickness of the quasi- 
2DEGs, and of phenomenological damping) in: (i) the RPA, which is expected to provide a 
good description when exchange and correlation are not important; (ii) the RPA with a K 
Hubbard local field correction to account for exchange- correlation effects (H-RPA); (iii) the 
self-consistent STLS approach, modified to include non-zero T e . Most previous theoretical 
work on exchange- correlation effects in electron bilayers has assumed T e = K (e.g. Refs. 
to ^), but a key feature of our experimental work is that we have raised T e substantially 
(~ Tp), and it is important to account for these observations. The capacity of the STLS 
approach to include exchange-correlation effects more reliably that the H-RPA at low T e is 
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not especially relevant here; indeed, at the relatively high 2DEG densities and momentum 
range considered here, the RPA and the STLS approach should produce similar results.! 
What is important is that the STLS can be readily modified to account for the effects of 
non-zero T e on the static local field corrections. 

As will be shown, there remain discrepancies between our experimental results and those 
predicted by these theoretical descriptions of electron fluids at non-zero T e , and this is one 
of the central conclusions of this work. In particular, the static STLS and more generally all 
static local field factors fail to describe effectively the effects of exchange and correlation in a 
non-degenerate 2DEG. However, our results do suggest that a dynamic local field correction 
could provide an effective description. 

The theoretical approaches employed are considered in Section |TJ and Appendix in 
which we also provide an overview of the RPA, H-RPA and STLS approximations for electron 
bilayers. The effects of non-zero temperature and 2DEG thickness on the STLS local field 
factors and pair-correlation functions are considered in Section [TTI|. The calculated forms of 
the lineshapes for Raman scattering from the bilayer plasmons at non-zero temperatures are 
discussed in Section [TV]. Details of the experimental measurements are provided in Section [V] 
and in Section [VI| we present the results of our Raman scattering measurements and discuss 
them in terms of the theoretical models. 

II. THEORY 

A. The RPA for an electron bilayer 

Our starting point for a theoretical description of the excitation spectrum of a double 
quantum well system is the RPA, which has been used extensively for the study of the collec- 
tive excitations in low-dimensional electron systems,&0 including electron bilayers.0S~l! 
Throughout, we will be concerned with excitations of angular frequency uo and wavevector 
(parallel to the planes of the 2DEGs) q. We assume that the two 2DEGs, which will be 
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labelled i and j, are isotropic. 

The charge- density fluctuation Spi(q,uj) in layer i induced by an external potential 
(f)j Cxt (q, to) in layer j is given by 

5pi(q,u) = Xij(q,u)(j>j ext (q,uj) , (1) 

where Xiji^i w ) are the elements of the bilayer density response function x(g, to). Within the 
RPA (and the STLS approximation), this is given by (see Ref. [|for further details) 
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where Vij(q) is an effective interaction between an electron in layer % and an electron in layer 
j; within the RPA this is the Coulomb interaction %(g). Xi(Qi u ) is the density response 
function for a single, isolated layer: 



x^r^Hj > ( 3 ) 

where \\ (<Z; u ) 1S the 2D Lindhard non- interacting electron polarisability for layer i. x°(<? ; w ) 
for a parabolic conduction band is given byil 

0, .... ^_2m /■ dk /°(k + q/2)-/°(k-q/2) 



X (q,u};n,T) = / — — ^ ? — . , (4) 

/r y (27r)^ k • q — m(u; + ij)/n 

where /°(k) is the Fermi-Dirac function for (temperature dependent) chemical potential /1 
and 7 is a phenomenolgical damping parameter which we take to be the inverse of the single 
particle relaxation time (7 = l/2r s ). Mermin has shown that for non- zero 7 the above form 
for x°(g,c<j) does not strictly conserve particle number;iil however, the sample quality in the 
present work was such that 7 was sufficiently small to ensure that the above approximation 
(Eq. fD is adequate. 

To evaluate x°(<Zj u ) f° r non- zero temperatures we follow the method set out by 
MaldagueH and rewrite the integral (Eq. |J) as:0 

f°° 1 

X (g,^/i,T)= / X \q,u-^T = Q) d// . (5) 

Jo 4fc B Tcosh KfjL - n')/(2k B T)\ 
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The polarisability at zero temperature with non-zero damping is given b; 



X °(q,u;;n = E F ,T = 0) 
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Using Eqs. ^ and 0, Eq. [| can be evaluated numerically to gives a form for the non- 
interacting electron polarisability at non-zero temperatures with damping. 

For a realistic system, it is also necessary to account for the finite spatial extent of the 
electron envelope wavefunctions along the confinement direction for the 2DEGs. This is 
achieved with the Coulomb form factors Fij, defined by Fij(q) = Vij(q)/v(q), where v(q) = 
e 2 /2eeoq is the 2D Coulomb interaction and e is the effective dielectric function for the media 
surrounding the bilayer. Interactions between plasmons and phonons were included using a 
frequency dependent dielectric constant.0 Using envelope wavefunctions ipi(z), determined 
(for sample structural parameters corresponding to the experimental sample discussed in 
Section |V|) from a self-consistent solution of the Poisson and Schrodinger equations (see Ref. 
21~1 for further details), the form factors 

Vij(q) 



v(q) 



dz / dz'e-^- 2 'l|<Mz)|V^')P 



(8) 



were evaluated numerically for wavevectors up to q = 6k F and fitted to functions of the 
forms 



F n (q) = l/(l + bq m : 
F 12 {q)=eM-d*q n ) 



(9) 



where b, d*, m and n are fitting parameters. We show in Fig. [I] the q dependence of Fn and 
F12 determined from Eq. || together with the fits obtained using Eq. M. The form factors 
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for the strictly two-dimensional approximation, F^(q) = 1 and F^(q) = exp(-qd), are also 
shown for comparison (d is the effective separation of the 2DEGs0); for large wavevector q 
the 2D approximations for the form factors fail significantly. 

The functional forms for the Coulomb form factors given in Eq. |9], together with 
X°(q,u>; H,T) from Eq. [5], provide x% (Eq- |3|) and %(g, u>; /i, T)(Eq. within the RPA. 

B. The optic and acoustic plasmons of an electron bilayer 

The energies of the collective charge-density excitations, the plasmons, of the electron 
bilayer system are given (in the absence of damping, 7 = 0) by the poles of the the density 
response function x(g, c<j). Direct insight into the physical character of the charge-density 
excitation spectrum is obtained by diagonalising x(g,cj), to give, as the diagonal elements, 
the density response functions for charge- density fluctuations in which the fluctuations in 
the two layers are in-phase (x+) and out-of-phase (%-) :@ 

x± = , 2 • (if) 

\xi X2J V V* 1 X2J 

The poles of x+(q,u) and X-(<?) u ) gi ye the optic plasmon (OP) and acoustic plasmon (AP) 
respectively. 

More experimentally realistic response functions, with non-zero damping (7 7^ 0) allow 
us to calculate the cross-section R(q, uj) for Raman scattering by charge density fluctuations 
from the imaginary part of the density response functions0 

R±(q,u) oc -(n(w) + l)Im[x±(g,w)] , (11) 

and hence the expected Raman lineshapes for the two plasmon modes; the energies of 
the plasmon modes are determined from the peaks in the intensities R±(q,u). The factor 
(n(u)+l) in Eq. [11], where n(u) is the Bose-Einstein occupation function, was not considered 
in Ref. ^D|, the result of which was to slightly underestimate the asymmetry and widths of 
the theoretically predicted Raman lineshapes. In this work we will not include this factor 
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in our theoretical analysis but rather correct experimental spectra for its contribution; in 
particular this allows for a more meaningful comparison between experiment and theory of 
peak widths, which give a measure of plasmon damping. 



C. The STLS approximation 

Although the RPA has been very successful in providing a description of the excitation 
spectra of 2DEG systems, it takes no account of the effects of exchange and correlation. This 
failing of the RPA becomes clear if we consider the instantaneous pair correlation functions 
gij(r), which give the probability of finding an electron at an in-plane position r in layer j, 
given that there is an electron at the origin in layer This is given by IS 

9 " (r) = 1 + 7^/ e "" |Sy(9) ^ yl (S? • (12) 

where iVj is the equilibrium number density in layer i. The static structure factors Sij(q) are 
related to the density response function u) by the fluctuation dissipation theorem:il 

oo 

SM = ^= V J coth (J?^j lm[ Xij ( q ,")} , (13) 

— oo 

where V denotes the principal value of the integral. 

The RPA gives the unphysical result of negative g%j{r) for small r, implying a much deeper 
exchange correlation hole than is physically realistic and hence overestimating significantly 
the correlation energy, the energy of interaction between the electron and its associated 
correlation holeE To correct for this, Singwi et al. allowed for a local depletion of the 
electron density around any given particle .EB For the two layer case this ansatz leads to an 
effective interaction between the responding electron in layer i and the induced charge in 
layer j of the form: 

Vij(q) =v ij (q) [1-Gij(q)] , (14) 
where the local field factors Gij(q) are given by 

G « w = vw/^Si [5y(|q " k|) "^' • (15) 
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For an electron bilayer system, the Cry(q) are generally positive and so the inclusion of 
exchange- correlation creates a 'hole' in the induced charge, reducing the strength of the 
interaction between the responding electron and the induced charge. 

From the local field factors Gij(q) we may determine u), the density response func- 
tion of the electron bilayer system, using the modified effective potentials Vij (Eq. [1^) in Eqs. 
H, [|and[l(| which in turn allows a determination of the structure factors Stj(q) through the 
fluctuation dissipation theorem (Eq. |13|). Thus the determination of the local field factors 
within the STLS approximation involves the self-consistent evaluation of Xij(<?> u ), <3ij(q) 
and G i:j (q) (Eqs. g and [ID . 

Our method for the self-consistent evaluation of these STLS local field factors Gij for a 
two layer system at finite temperature is described in detail in Appendix and mirrors the 
approaches of Tanaka and Ichimaru, who studied a one component plasma system at finite 
temperature in three dimensions within the STLS, principally with a view to modelling stel- 
lar cores,ii and of Schweng and Bohm who extended this work to examine a one component 
electron liquid in two dimensions.il Note that for simplicity we assume that the two layers 
are of equal density i.e. Ni = Nj = N. This approximation is justified for the sample studied 
experimentally, where any asymmetries between the wells were minimal; calculated Fermi 
energies for the two wells were within 0.5% and no experimental asymmetries were observed 
in, e.g. photoluminescence measurements or in the dispersive behaviour of the acoustic and 
optic plasmon modes at low temperature. 

The local field factors evaluated within the self-consistent scheme outlined in Ap- 
pendix [A] are then used to determine the diagonalised density response functions x±( < li UJ ) 
(Eq. pi]), determined with the Lindhard polarisability given by Eqs. |5] and |6[. The RPA 
result is recovered if the local-field factors Gij(q) are set to zero. These density response 
functions allow us to determine (through Eq. [11]), as a function of electron temperature T e , 
the Raman scattering lineshapes for the acoustic and optic plasmons of the electron bilayer, 
which we then compare with our experimental results. 
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D. The Hubbard approximation 



H-RPA 



The simplest way to go beyond the RPA is to avoid the self-consistent scheme set out 
in the previous Section and to use instead the Hartree-Fock approximation for the static 
intra-layer structure factor Su(q)p which accounts only for the presence of the Pauli hole 
around each electron at T = 0. The effects of non-zero temperature and contributions to 
the local field from intra- and inter-layer Coulomb correlations are therefore neglected in 
this approximation — the Hubbard correction. The resultant integral for Gij(q) (Eq. |T3|) 
can then be evaluated explicitly (within an approximation) to yield:0 

Gf 3 (q) = J 2 vV + Wi 1 1 3 (16) 
{ ifz^j 

where k-Fi is the Fermi wavevector in layer i. 

This Hubbard local field correction G^(q) was used to account for the effects of electron 
correlations in the analysis of the Coulomb drag measurements of Ref. and in our pre- 
liminary report on the temperature dependent Landau damping of the acoustic plasmon in 
an electron bilayer.0 



III. CALCULATIONS OF CORRELATIONS IN ELECTRON BILAYERS 

We present in this section the dependence of the STLS local field corrections Gy(q) and 
pair correlation functions g%j{r) of an electron bilayer on temperature and 2DEG thickness. 
For all the calculations presented here, we have taken the parameters for the electron bilayer 
to be those for the sample investigated experimentally (see Section |V|). In fact the results 
show that the effects of temperature and non-zero layer thickness on the inter- and intra-layer 
STLS local field corrections and correlations are not significant for this sample. 
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A. Temperature dependence of local field factors 



Our calculations of the temperature dependent local field factors Gij(q) for electron 
bilayers reproduce the results of earlier calculations in the limit of zero temperature and 
layer thickness. In particular we see similar behaviour of the local field factors as functions 
of number density N (or equivalently r s = (iWal)" 1 / 2 ) and interlayer separation (d), to 
that described by Liu et aM 

In Fig. |2] we show the q dependence of the local field corrections Gn and G12 as a function 
of temperature {6 = T e /T F ). Although the wavevector range (to 6 hp) is not relevant to our 
experimental work, it was necessary to determine Gij for large q in the self-consistent STLS 
scheme (Eq. [15]), and we present these results for completeness and to illustrate the integrity 
of our calculations. The dependence on temperature of Gn and G 12 for q = 0.2 k F (the 
typical magnitude of wavevector accessible in the Raman experiments) is illustrated in Fig. 
[3]. Although G12 increases significantly with temperature, the magnitude of this inter-layer 
local field factor is negligible for this wavevector. Gn, on the other hand, is much larger 
and is reduced (although by only ~ 5%) on increasing the temperature from absolute zero 
to the Fermi temperature. 

Fig. |] shows the dependences of gu(r) and guir), the pair correlation functions, on r, for 
several temperatures; the minimum in the pair correlation functions describes the exchange- 
correlation hole. The dependence on temperature of the intra-layer pair correlation function 
g\i{r) is similar to that for a single layer JH as T e increases the radius of the exchange- 
correlation hole is reduced, accompanied by a slight deepening of the hole, in accordance 
with screening sum rulesJH 



The inter-layer pair correlation function gu{r) for these comparatively high densities (r„ ~ 
1.3) is, as expected, close to 1. 




(17) 
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B. Quantum size effects 



In Fig. |5] we show the q dependence of the intra- and inter-layer local field corrections 
G\i an d G12 determined with the finite thickness effects of the quantum wells taken into 
account (using the functional forms for Fu and F12 given in Eq. ||), which we compare with 
the corresponding results determined within the strictly two-dimensional approximation (for 
which F^q) = 1 and F® 2 {q) = exp(-qd)). Our results are in line with those for one layer 
systems,^ in that correlations are reduced in going from the ideal to the quasi-2D regime. 
The greatest discrepancies occur for large wavevector q, where the 2D approximations for the 
Coulomb potentials fail significantly (see Fig. [l]). In particular, the inter-layer correlations 
for a real system do not continue to increase as a function of q, but instead reach a maximum 
and then decline. 

Fig. H shows a similar comparison between the quasi-2D and ideal-2D cases for the 
intra- and inter-layer pair correlation functions, gn and gi2- Compared with the ideal-2D 
case, electron-electron effects are expected to be reduced in the quasi-2D regime because 
the spread of charge in the ^-direction reduces the intra-layer Coulomb interactions; so 
there should be a shallower intra-layer exchange correlation hole in the quasi-2D regime, as 
confirmed in the calculations. This shallowing of the intra-layer hole is accompanied by a 
deepening of the inter-layer hole, as electrons in different layers do not in general want to sit 
opposite each other. These effects of the charge spreading in the ^-direction are significant 
for this density (r s ~ 1.3) because the mean distance between the electrons (~ 1/kp ~ 90 A) 
is significantly less than the layer width. 

IV. RAMAN SCATTERING FROM CHARGE-DENSITY FLUCTUATIONS 

Fig. [7] shows calculated spectra for Raman scattering by the AP (derived from 
—lm[x-(q,w)]), the OP (from —Im\x+{q,w)]) and the SPE (from — Im[x°(g, &)]) for T e = 10 
K and 50 K. As T e increases, the SPE spectrum spreads upwards in energy, as expected. 
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This results in an increasing population of electrons travelling with the same velocity as the 
phase velocity of the AP. These electrons are able to exchange energy efficiently with the 
plasmon and so enhanced scattering of the AP by SPE occurs, resulting in an asymmetric 
broadening of the AP peak on the low-energy side, a clear signature of Landau damping.0 
The degree of AP broadening, and its asymmetry, can be characterised (this will be dis- 
cussed later) in terms of the half-widths on the higher and lower energy sides of the peak. 
In these calculations the occupation factor (n(o>) + 1) in the Raman intensity (Eq. [ll]) has 
been omitted, to highlight the effects of Landau damping. The effect of this factor, which 
becomes more significant with increasing temperature, is to increase the Raman intensity at 
low Raman shifts, leading to a further slight increase in asymmetry of the AP peak. 

Note that the OP peak, at rather higher energy, does not overlap the SPE continuum, 
and is not detectably broadened by Landau damping, even at 50 K. Fig. |7| also shows the 
upward drift in the energies of the AP and the OP with increasing temperature, a well 
known classical effect to which we will return. 

V. EXPERIMENTAL MEASUREMENTS 
A. Sample details 

The sample investigated (T229; sample C in Refs. |2D] and [H]) consists of two GaAs 
quantum wells of width L w = 180 A, separated by an Alo.67Gao.33As barrier of width L h = 
125 A, giving an effective inter-2DEG separation d ~ 305 A. These parameters are close 
to optimal for this study, in that Lb is large enough to preclude inter-well tunnelling^ (we 
have throughout assumed a quantum mechanically decoupled system), while the inter- layer 
separation d is small enough to lead to significant inter-layer Coulomb interactions, which 
depress the energy of the AP mode towards the SPE continuum and increase its susceptibility 
to many-body effects. 

The number densities N\ and N2 in the two wells under laser illumination were obtained 
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from a fit, calculated within the STLS approximation, of the low-temperature experimentally 
determined AP and OP dispersions. HI It was found that Ni and N 2 were equal to within 
~ 5%, as expected since the two wells were grown to be equivalent. N\ — N 2 — (1.95 ± 
0.10) x 10 11 cm~ 2 ; see Ref. ^1] for further details. 

The sample was mounted in a He-cooled cryostat and its temperature controlled by 
electrical heating. 



B. Determination of the electron temperature 

The effect of the 2DEG temperature T e on the plasmons was studied by varying Tl, 
the lattice temperature. An accurate determination of T c was critical for the interpretation 
of our experiment, and was achieved by carefully studying the quantum well and bulk 
photoluminesence (PL) lineshapes, determined under the same conditions as the Raman 
scattering measurements presented below; i.e. for excitation photon energies > 1.65 eV and 
high powers (40 W cm -2 ), with T e > 25 K. The strength of the bulk PL signal from the 
GaAs buffer layer was found to be some 25 times greater than the quantum well signal, 
as shown in Fig. |8](a); the lower energy signal is due to the bulk GaAs, and the weak 
shoulder at ~ 1.53 eV is the quantum well PL (this assignment has been confirmed by PL 
and PLE measurements performed with lower power densities and electron temperatures). 
Consequently the PL tail from the bulk GaAs must be fitted and subtracted from the 
experimental spectra to give the quantum well PL; a fit to the bulk PL lineshape is shown 
in Fig ||(a). 

Having stripped the PL spectrum of the bulk PL, the remaining quantum well PL line- 
shape £qw(^) was taken to be: 

C QW {u, T e ) oc /eCEe,k, T e )f h (E hjk , T e )5{E' g + £ c , k + £ h , k - Hu) , (18) 

k 

where / e (h) are the Fermi-Dirac occupation factors for electrons (holes) with kinetic energy 
E'e.k (-E^k) and wavevector k. E' is the quantum well band-gap and ftju is the PL photon 
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energy. Given the temperature range over which the measurements were made, simplifica- 
tions to the form of £qw(w,T e ) (e.g. by making assumptions about the degeneracies of the 
electrons and holes) were not possible, and the experimental quantum well PL lineshapes 
were fitted to £qw(w) assuming only that: (i) Coulomb scattering renders the electron and 
hole temperatures equal; (ii) the total extrinsic electron density (1.95 x 10 11 cm -2 for each 
quantum well) did not change with temperature; (iii) the photocreated and thermal elec- 
tron and hole population densities (used as fitting parameters) were equal. The electron 
and hole effective masses were taken to be 0.07 and 0.18 respectively. We show in Fig. ||(b) 
a logarithmic plot of the experimental QW PL for three different temperatures, along with 
fits of the form given in Eq. [18|. The data does not fit a straight line (the exponential form 
often used) except well into the higher energy part of the PL tails where the errors in the 
data are significant. 

By fitting the bulk and QW PL lineshapes in this way, we were able to determine the 
temperature T e of the 2DEG for each Raman measurement with a maximium error (for the 
highest temperatures) of ~ 5%. We were able to derive a value for the lattice temperature 
Tl from the form of the variation of E g (determined from the bulk PL) with Tl:0 

F-F(T-H) (5-408 xlO- 4 eV)(T L /K) 2 

E g - E g (T L - 0) (Tl/K ) + 204 (19) 

We show in the inset of Fig. |](a) the variation of T c with the lattice temperature Tl; the 
laser illumination is seen to heat the electron gas substantially above the lattice temperature 
for low Tl (due to the relatively high photon energy and power densities used), with T e and 
Tl converging for Tl > 30 K. 



C. Raman scattering measurements 

Angle resolved Raman scattering measurements were made under incoming resonance 
conditions, and with the polarisations of the incident and scattered light parallel to allow the 
observation of plasmons.0 The band-gap decreased with increasing lattice temperature Tl, 
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so to maintain the conditions for an incoming resonance the laser wavelength was increased 
correspondingly. For the low-temperature measurements the energy of the resonance used 
was 1.656 eV. Rather than determine a resonance profile for each measurement temperature, 
the shift with increasing T e in the resonance energy from this value was determined using 
Tl and the known variation of the semiconductor band-gap with temperature (Eq. [T9| ). 

As the temperature increases the tail from the bulk bandgap PL increases considerably, 
and can overwhelm the comparatively weak Raman signals. Fortunately the resonance 
utilised here was sufficiently far from the bandgap PL to allow Raman measurements to 
be made over a reasonably large temperature range. The band-gap PL, and the hot PL 
originating from the inter-band transition responsible for the Raman resonance, limited the 
Raman measurements to T e < 110 K. 



D. Single particle lifetime 

In our calculation of the Raman scattering lineshapes for the plasmon modes, the phe- 
nomenological damping constant 7 (see Eq. |]) is taken to be the inverse of the 2DEG single 
particle relaxation time r s , which we determined experimentally from the dependence of the 
magnitude of Shubnikov-de Haas oscillations with magnetic field; only one Shubnikov-de 
Haas oscillation frequency was observed, indicating that the two quantum well densities 
were the same. The amplitude (Ai?) of Shubnikov-de Haas oscillations is given by:! 

AR = 4R X(T) exp(- — ) , (20) 

where Rq is the zero field resistance, uj c the cyclotron frequency, and X(T) a thermal damping 
factor given by: 

= 2n 2 k B T/hw c 
1 ' sinh(2vr 2 A;BT/^ c ) ' 1 ' 

For our samples, r s was found to be of the order of 3.5 ps for Shubnikov-de Haas oscillations 

measured at 4 K. Using this value, we were able to fix the low temperature half width of 

the plasmons as frry = fr/2r s = 0.09 meV. 



VI. RESULTS 



A. Temperature dependent Raman scattering 

A Raman spectrum for T e = 93 K is shown in Fig. 0(a); at such high electron temper- 
atures the Raman signals are superimposed on quite intense background PL arising from 
the quantum well inter-band transition responsible for the Raman resonance ('hot' PL, at 
small Raman shifts) and also from the GaAs buffer layer (bulk PL, at high Raman shifts). 
In addition, the importance at such elevated temperatures of the occupation factor in the 
Raman intensity (n(u) + 1 in Eq. [TT| ) is also responsible in part for the large intensity at 
small Raman shifts. In order to facilitate comparison between experiment and theory, espe- 
cially for the AP damping as characterised by Raman peak widths, the experimental spectra 
I{oj) were first divided by the occupation factor (n(u) + 1), as shown in Fig. (|(a). Fits to 
the lineshapes I(uj)/(n(uS) + 1) were then obtained using exponential tails to describe the 
PL contributions to the experimental lineshapes; the result of such a fit is indicated in Fig. 
0(a). The exponential bulk and hot PL tails were then removed to give a corrected Raman 
signal, also shown in Fig. |9|(a), in which the effect of the occupation factor (n(uj) + 1) has 
been removed. 

Fig. 0(b) shows Raman spectra measured for a range of electron temperatures T e , for 
an in-plane wavevector transfer q = 1.6 x 10 5 cm -1 . These spectra have been corrected as 
described above to leave signals due principally to the plasmons. Two peaks are present, 
which we assign to the acoustic and optic plasmons of the 2DEG@ The OP is seen to 
broaden symmetrically as T e rises, whereas the AP is seen to broaden asymmetrically and 
to a much greater extent. 

We ascribe the asymmetric broadening of the AP to Landau damping arising from the 
increasing interaction with the lower lying SPE continuum as T e is increased (as discussed 



in Section IV and illustrated in Fig. m). The broadening with T e of the OP peak may 



be attributable to increased scattering by unionised impurities in the sample bulk as the 
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temperature increases; this scattering is greater for the OP than the AP as away from the 
bilayer region the electric fields associated with the symmetric OP are much greater than 
those of the antisymmetric AP.il We have also shown that this effect can account for the 
observed dependence on 2DEG mobility of the relative strengths of Raman scattering by 
the AP and OP.0 

The observed AP width involves several contributions — Gaussian terms such as inhomo- 
geneous broadening effects, spectrometer resolution, and broadening due to the spread of in- 
plane wavevectors accepted by the numerical aperture of the collecting lens, and Lorentzian 
contributions from electron scattering (described by the parameter 7 in Eq. ^) and Landau 
damping. These combine according to: 

Au L Au tot + (Au G ) 2 = Au 2 tot , (22) 

where Au^Auq) is the Lorentzian (Gaussian) contribution to the total half- width at half 
maximum (HWHM), Aa> tot @ To allow easy comparisons with the calculated spectra, the 
Gaussian contributions to the half-widths were deconvolved from the experimental spectra 
by adapting a method originally due to Dobryakov and Lebedev.il As a result the net 
Gaussian contribution to the FWHM (full-width at half maximum) was estimated to be 
0.06 meV, in good agreement with the observed width of the Rayleigh scattered laser line 
(cf. the single particle relaxation FWHM 27 = 0.18 meV). 

Because the AP lineshape is asymmetric, HWHMs on either side were measured and Eq. 



22] used to subtract Aujq in each case to produce yi and y 2 , the HWHMs of the Lorentzians 
on the low and high energy sides, respectively (illustrated in Fig. |^(a)). The width of the AP, 
once Gaussian contributions have been subtracted, is a measure of the total homogeneous 
damping experienced by the plasmon mode. The difference between yi and y 2 , on the other 
hand, gives information about the energy-dependent damping responsible for the observed 
asymmetry of the AP peak; thus y\ — y 2 gives a direct measure of the Landau damping 
experienced by the AP. In this way, by comparing y± and y 2 with theoretically calculated 
half-widths, we have a direct handle on the many-body interactions of the AP with the SPE 
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continuum. 

We should note that the SPE band itself is potentially observable in both polarised and 
depolarised Raman spectra. However, at low temperatures no SPE signal was observed in 
either polarisation for q = 1.6 x 10 5 cm -1 . For q = 1.16 x 10 5 cm -1 the SPE band is in a 
spectral region too close (a Raman shift less than ~ 1 meV) to the exciting laser line to be 
observed. 



B. Total AP damping 

We now consider the comparison between the experimentally observed and theoretically 
calculated AP total linewidths as a function of temperature. Fig. |K](a) shows a plot of 
the experimentally determined FWHM (yi + y 2 ) for q = 1.16 x 10 5 cm -1 , after removing 
the Gaussian contributions as described above. The experimental uncertainties in these 
FWHM values are of the order of ± 0.02 meV for low T e , rising to ± 0.06 meV for higher 
temperatures. Fig. 0(b) shows the corresponding results for q — 1.6 x 10 5 cm -1 . The three 
curves show the results of calculations of the FWHM determined within: (i) the RPA; (ii) 
the H-RPA; (hi) the finite temperature STLS formalism described in Section ||. 



We can see from Fig. |10| that the larger the magnitude of the local field correction G n , 
the greater the increase in damping with increasing T e ; as the Hubbard intra-layer local field 
correction G n is greater than that obtained within the STLS, for a given T e the H-RPA 
predicts a greater width than that calculated using the STLS approximation, which is in 
turn greater than that obtained within the RPA. 

For low T e , for both wavevectors, the experimentally determined FWHM and the calcu- 
lated values agree quite well, especially as it should be noted that no free fit parameters were 
used for these curves — the K theoretical FWHM was fixed from the experimental single 
particle relaxation times (Section |V D| ). But since the experimental uncertainties at these 
values of T e are smaller than the point size in the Figure, the better agreement between the 
data and the H-RPA and STLS curves, as opposed to the RPA calculation, is significant, and 
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suggests that it is necessary to include exchange-correlation effects even for these densities. 

However, the experimentally determined plasmon widths increase with T e much more 
slowly than is predicted theoretically, and for T c > Tp drop below even the values predicted 
by the RPA. This could be interpreted in terms of a reduction in the significance of correlation 
effects at these temperatures, at which the RPA then provides a better description. However, 
the STLS calculation, which should take full account of the effects of non-zero temperature, 
does not show this behaviour, reflecting the relative constancy of the intra-layer local field 
correction Gn over this temperature range (see Fig. ||). So it appears that for T e > Tp, 
further theoretical developments are necessary. 

These calculations take no account of any dependence on temperature of the single 
particle relaxation rate 7. However, if anything we would expect the single particle lifetime 
to decrease with increasing temperature, leading to even larger theoretical linewidths than 
those calculated here and increasing yet further the discrepancy between experiment and 
theory. 

C. AP asymmetry 

For a given T e it is expected that the low energy side of the AP peak will be more affected 
by Landau damping than the high energy side, as the SPE continuum lies below the AP 
dispersion curve. A measure of this asymmetric broadening is (yi — 1/2), which is plotted in 
Fig. [TT] for both experiment and theory, again for both q = 1.16 and 1.6 x 10 5 cm -1 . 

It is expected that (2/1 — 2/2) —> as T e — > (the peak shape would be a symmetric 
Lorentzian) and increase with increasing T e . Experimentally, for q = 1.6 X 10 5 cm -1 (Fig. 
|TT| (b)), (yi —2/2) increases with low T e as expected, but at higher temperatures (T c ~ 70 K) 
it becomes roughly constant. For the smaller wavevector of q = 1.16 x 10 5 cm -1 (Fig. [□] 
(a)), the discrepancy between experiment and theory is greater and (yi —2/2) is essentially 
zero over the measured temperature range. In contrast to this, all the theoretical models 
predict a continuous increase in {yi — y 2 ) with rising T e . It is important to note that any 
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temperature dependence of the single particle lifetime (~ I/27) is unlikely to affect this 
measure of the asymmetry of the peak; we are sensitive here only to the effects of Landau 
damping. 

Another measure of the Landau damping is provided by the degree of asymmetry, a = 
(Vi — 2/2)/ (2/1 + 1/2)1 shown in Fig. |T2| for experiment and for our STLS calculation. For 
the two wavectors considered, the STLS theory predicts essentially the same variation with 
temperature, with a levelling off near Tp. This insensitivity of a to the wavevector is 
to be expected, as both the AP and the SPE continuum have a similar, approximately 
linear, dependence on q, for the small wavevectors (q -C /cf) considered here. However, the 
experimental results for the two wavevectors are markedly different. For q = 1.6 X 10 5 cm -1 
there is a fairly good agreement between experiment and theory over the whole temperature 
range, although the experimental uncertainties are large (~ ± 40% for high T e ); note, though, 
that the calculations overestimate the FWHM of the AP peak at higher temperatures. For 
q = 1.16 x 10 5 cm -1 there is a significant discrepancy between experiment and theory. 



D. AP energy 

Figs. [13| (a) and (b) show the energy c^ap of the acoustic plasmon as a function of T e 
for q = 1.16 and 1.6 x 10 5 cm -1 , for both experiment and theory. The 2DEG density for 
the calculations was set so that the AP dispersion obtained within the STLS approximation 
matched the experimentally determined dispersion at low temperature.0 (Previously^ the 
density was obtained using a fit determined within the H-RPA, resulting in a slightly higher 
density than that used here; this results in the slight shift of the RPA and H-RPA curves 
with respect to our data for (yi +1/2) (Fig. [H]) and (yi — y%) (Fig. [□]), compared with those 
presented in Fig. 4 of Ref. |2"0| .) 

All the theoretical curves show a large increase (~ 30%) of the AP energy with tempera- 
ture, for both wavevectors. The OP energy c^op was also found to increase with T e , although 
by only ~ 10%. Such behaviour is well known for classical Maxwell-Boltzmann plasmas,@ 
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and similar calculations carried out for a 2D plasma correspond closely to our RPA results 
at higher temperatures. 

Fig. [13] shows that the inclusion of exchange-correlation effects reduces the calculated 
values of uap, as would be expected from the softening of the effective interaction potentials 
Vij (Eq. [TJ). The weak temperature dependence of the STLS local field correction (see Fig. 
|H) does not significantly alter the variation of the plasmon energy with T e . 

Experimentally, for T e < 50 K (Tf = 78 K), uap increases with T e , in keeping with the 
calculations. However, for both wavevectors, for T e > 50 K, ujap remains approximately 
constant, in contrast with the theories, all of which predict a steady increase in the plasmon 
energy. The variation of the OP energy cuqp with temperature follows essentially the same 
pattern as that of the AP (see Fig. |9|(b)), although there is better agreement between 
experiment and theory due to the smaller predicted increase of u>op (~ 10%) with T e . 

E. Discussion 

Above T e ~ 50 K, there are significant discrepancies between experiment and theory. All 
three calculations show an increase in uj^p, contrary to the observed plateau behaviour, and 
all three overestimate the width and asymmetry. Not only does the RPA fail to describe our 
results, but so does the full temperature dependent STLS scheme. We believe that these 
discrepancies are significant, suggesting a failure of the static local field theories used here to 
describe adequately the excitation spectra of non-degenerate low- dimensional electron gases. 

To demonstrate that this is the case, we must firstly to examine whether alternative 
causes of the observed plasmon behaviour, such as changes with T e of the 2DEG densities, 
or of their background screening, could explain our experimental results. Fig. [14] explores 
the effect of possible changes in the 2DEG densities by plotting, for various 2DEG densities, 
calculations of c^op, <^ap and the AP FWHM Suap = (yi + Hi) as functions of T e for 
q — 1.6 x 10 5 cm -1 , together with the corresponding experimental results. (For expediency, 
the STLS local field factors self-consistently determined for a density of 1.95 x 10 11 cm -2 (Fig. 
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were used in the calculations for all densities; this can be justified by the small variation in 
the low T e values for Gn and G 12 over the relatively small range of densities considered.) The 
observed temperature dependence of uj ap (Fig- |H|(b)) would require a significant reduction 
in 2DEG density (by about 30% for the highest measurement temperatures), though a much 
smaller reduction (~ 15%) would account for the behaviour of c^op (Fig. |TJ](a)). Moreover, 
Fig. |l4l(c) shows that the AP damping 5u>a_p actually increases with decreasing density, 
worsening the discrepancy between experiment and theory; for a given T e , the smaller 
resulting from a decrease in density effectively increases the reduced temperature T e /Tp and 
enhances the Landau damping. So although a reduction in density with temperature could 
account for the observed T e -dependence of ujap, the corresponding temperature dependence 
of co>op an d Sujap cannot be reconciled to this. 

We have assumed throughout that the two 2DEGs are of equal density, an assumption 
borne out at low temperatures from the form of the plasmon dispersionsBEll However, a 
temperature induced asymmetry in the electron distribution between the two quantum wells 
(the total electron density remaining constant) might occur; calculations (within the H-RPA 
for computational simplicity) show (Fig. |T5| (b)) that this would result in a reduction of ujap, 
and that agreement between experiment and theory is possible at the highest measurement 
temperatures for N1/N2 ~ 3 (i.e. N\ ~ 2.9 and N2 ~ 1.0 x 10 11 cm -2 ). However, not 
only is this a very significant asymmetry, but again neither the T e -dependence of u>op, or 
that of Suap, fall into line with this. An asymmetric density distribution results in a slight 
increase in cuqp (Fig. [15] (a)) as the two plasmon modes decouple and the OP becomes more 
localised in the higher density layer; at the same time the AP becomes more localised in 
the lower density layer, leading to the observed reduction in energy. An increasing density 
asymmetry also leads to an increase in the Landau damping of the AP, which comes closer 
to (and enters, for a large difference between the densities) the SPE continuum of the higher 
density layer.i 

The variations of uqp and uap with T e (see Figs. [TJ] and [T5|) also indicate that changes 



in background screening (from, e.g., un-ionised donors in the doped AlGaAs layers) are not 
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responsible for the reduction in uap from the theoretical values. If this were the case, we 
would expect ujqp to be affected to a much greater extent, as seen in the widths of the two 
plasmon peaks (see Section |V1 A| and Fig. ^). In addition, although increased screening 



would reduce u\p, the SPE continuum would be unaffected, resulting in increased Landau 
damping, contrary to observation. 

Therefore it appears that the inability of our calculations to model our experimental 
results for the energy and damping of the AP cannot be simply attributed to an experimental 
artifact, such as a decrease in 2DEG density with temperature. We must thus examine 
whether the use of a static local field correction provides an accurate description of the 
effect of exchange and correlation on the excitation spectra of non-degenerate quasi-two- 
dimensional electron systems. In order to obtain better agreement with experiment for 
the AP energy, the magnitude of Gu must increase with temperature, whereas to obtain an 
improved fit for the measured AP damping, Gn must decrease.0 This contradiction suggests 
that a real local field factor, (as calculated within the framework of a static theory) cannot 
explain our observed results. This is one of the central conclusions of this work. 

We are led therefore to consider whether a complex local field correction, Gn = G' X1 + 
iG'd, as suggested in the context of a dynamic local field theorySil is any more successful 
in characterising the observed behaviour. For a given temperature T e , G' u (the real part 
of Gn) was determined from the energy of the acoustic plasmon c^ap; G'u (the imaginary 
part of Gn) was found to have a negligible affect on uj^p- Singwi et a/.i have argued that 
the imaginary part of the local field correction is important in determining the lifetimes 
of plasmon modes. So, using the derived value for G' xl , the imaginary part G' n was then 
determined by fitting the plasmon damping 5a; ap- In order to give a reduction in the AP 
damping, to bring theory in line with experiment, G' n was determined to be positive. As well 
as now providing complete agreement between experiment and theory for the T e -dependence 
of cuap and 8ui ap, the resulting complex local field factor, Gn = G'n + iG'd, was also found to 
provide a more accurate description of the optic plasmon energy uop and the AP asymmetry 
for both wavevectors q considered here. The variation with temperature T e of these fitting 
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parameters, the real and imaginary parts of Gn, is shown in Fig. [16| for the two wavevectors 
q. For T c < Tp/2, the real and imaginary parts of the local field correction are similar to 
the values obtained within the STLS; in particular G'[ x is zero. However, for T e > Tp/2, the 
temperature range over which the static STLS approach was found to fail, we see that both 
the real and imaginary parts of the intra- layer local fields increase with T e . This is perhaps 
surprising, given our initial expectation that exchange-correlation effects should become less 
important with increasing temperature. Nevertheless, it is significant that the only way in 
which we have been able to obtain agreement between experiment and theory is through 
the incorporation of a complex local field correction. 

Finally, we observe that our conclusions are consistent with the results of measurements 
of the AP-mediated Coulomb drag between two electron layers in samples similar to that 
studied here. EE In particular, in such experiments the onset temperature for plasmon 
enhancement of the drag rate is determined by the energies of the acoustic plasmons. A 
suppression of the AP energies from that expected within the RPA, as observed in our 
Raman measurements, should lead to the plasmon enhancement in the Coulomb drag rate 
occurring at a lower temperature than that predicted by the RPA; this is indeed observed 
in Refs. 123 and ESI. On the other hand, an increase in the damping of the plasmons leads to 



an increase in the magnitude of the plasmon enhancement in the drag rate. Hill et al. found 
that although the use of a Hubbard static local field factor resulted in better agreement 
between theory and experiment for the onset temperature, due to the suppression of the 
AP energies, the resultant increase in the plasmon damping leads to an overestimate of the 
magnitude of the enhancement ,@ which mirrors the results of the Raman measurements of 
the AP damping presented here {i.e. at higher temperatures we find the plasmon damping is 
reduced from that predicted by theory). Our use of a complex local field factor to reconcile 
the behaviour of the plasmon energy and damping is consistent with the conclusions of Noh 
et al, who suggested that a dynamic local field theory could account for their Coulomb drag 
results,il although they implied that this would lead to a broadening of the plasmon modes, 
whereas we find the damping (as determined from the Raman lineshape) to be less than 
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that predicted by the RPA. 



VII. SUMMARY 

To conclude, we have analysed the temperature dependence of the Landau damping of the 
acoustic plasmon in an electron bilayer, with a view to obtaining a measure of the importance 
of exchange-correlation effects in such systems. To this end, we have modelled the density 
response of our system using the Singwi-Tosi-Land-Sjolander (STLS) approximation, in 
which we have taken full account of the effects of finite temperature. The most significant 
message of these calculations is that non-zero temperature correlations within the STLS 
differ little from their zero temperature counterparts at low wavevector. Increasing the 
temperature leads to an increase in the kinetic energy of the electron gas, implying a decrease 
in the significance of many body effects in determining its behaviour. Naively we could have 
expected our temperature dependent local field corrections to show a transition from the 
H-RPA (a model including exchange-correlation effects to some degree) to the RPA (a mode 
in which exchange and correlation are not included) .0 If this were the case we would have 
expected a much greater decline in the local field corrections towards zero as T e is increased. 
However, this does not seem to be the case for the static local field corrections we have 
calculated within the STLS. 

Experimentally, we determined the electron temperature of the 2DEGs from measure- 
ments of the photoluminescence, which gives an accurate determination of the temperature 
under the same conditions as the Raman measurements. Our Raman lineshapes displayed 
the clear signature of Landau damping of the acoustic plasmon, indicated by the asymmet- 
ric form of the Raman peak. Careful analysis of the spectra allowed us to extract three 
parameters which characterise the lineshape and energy of the acoustic plasmon: y% and 
U2, the low and high Raman shift half widths of the plasmons, and the AP energy itself. 
These parameters were compared with theoretically calculated values, determined within 
the RPA and within the H-RPA and STLS local field theories, thereby giving us a means to 
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analyse the success of these theories in modelling electron gases of intermediate degeneracy. 
We found that for T e < Tp/2 the inclusion of exchange- correlation effects using a local field 
correction is necessary — indeed sufficient — to obtain good quantitative agreement for 
the asymmetric Landau damping. Conversely, for T e > Tp/2 none of the theories modelled 
successfully the AP asymmetry, or its upward shift in energy with increasing T e . Further, 
the behaviour of the AP energy and damping cannot be reconciled to the calculations by 
a reduction with T e of the 2DEG density, or a redistribution of electrons between the two 
quantum wells. This leaves only exchange-correlation effects as the cause of the discrepancy 
between experiment and theory. We find that the STLS approach and indeed any static 
local field theory cannot effectively model the exchange-correlation effects apparent in the 
experimental data. 

We should note also that we observed a greater discrepancy between experiment and the- 
ory for a wavevector of q ~ O.lfcp than for q ~ 0.14/cp. Indeed, for the lower wavevector the 
Landau damping of the AP appeared to be effectively constant over the whole temperature 
range considered. 

There are a number of alternative, although less widely applied, theoretical approaches 
to that of the STLS. It may be that a non-local scheme (see, e.g., Refs. [TT] and |T^) may 
be necessary for a true description of the excitation spectra of low-dimensional electron 
systems. However, from our work the most promising candidate is probably a quantum 
STLS approach, employing a frequency dependent local field correction,0'@ which may go 
towards providing a more complete agreement between experiment and theory; we have 
shown that we are able to obtain excellent agreement between experiment and theory using 
as a fitting parameter a complex local field factor, as suggested in the context of a dynamical 
local field theory. Our work shows that the temperature dependence of Landau damping is 
likely to be a valuable experimental probe of these theories. 
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APPENDIX A: 

We describe here our implementation of the STLS scheme for electron bilayers. The 
density response functions Xij(li u ) w iU have poles on (or infinitesimally close to) the real 
uj axis which correspond physically to the energies of the plasmon modes, while the poles of 



the coth function in Eq. [13] will lie equispaced along the imaginary u axis. A Wick rotation 
(u — > uS) allows us to evaluate the integral for the static structure factors Sij (q) in Eq. [13] as 
a contour integral along the imaginary uj axis, where the poles of the integrand are now just 
the poles (j v = 2nii>k^T/h, v = 0, ±1, ±2...) of coth(^o;/2A;BT). The residue theorem then 
allows a transformation of the integral into a sum over imaginary (Matsubara) frequencies 
giving: 

u=oo 

Sij(q) = — y== Yl XiAvnu) ■ (Ai) 

For convenience, we now assume that the two layers are of equal density, i.e. Aj = Nj = 
N. Defining $y(g,7„) = -(E F /N)xij(q,lu), we have: 

u=oo 

SiM= E ^iu)e , (A2) 

where 6 = k-^T / E-p = T/Tp and Ep and Tp are the Fermi energy and temperature. From 
Eq. S we have: 



i - vys* 



( 1 ~ * N 

1 -«i2$l 
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(A3) 



<&2D has been evaluated explicitly by Phatisena et al. who derived the form:0 

•a.h.7.) - [ *m [{ql/i _£™ { ^ +q ^ v/1 « , ( A5) 

where /(/c) is the Fermi-Dirac distribution function and 

t-(20) = qij ,%^ ■ (A6) 



Note that in the case of a negative denominator 20 has to be replaced by n — 2<p. 

Thus in order to evaluate Sij(q) we must evaluate Eq. |A2] using Eq. [A~3| . 511(g) can be 
written as: 

Sn® = £ : , , (A7) 

which is reformulated as 

e f\ a <h ( \ , ^^2(9)*i(9,7.) / As x 

S11 (?) = V ^iW,7. + V 1 ~2 / n ^2/ V • A8 

^ ^ 1 - v f, <Pf a, 7„ 



f=— OO 1/= — 00 



S 2D (g) AS 
The first term in Eq. [A8| is the structure factor for a single electron layer, S^d^)- The 
second, interlayer term AS is 0($f) and acts as a (small) perturbation, reflecting the small 
extent to which the second layer affects intra-layer correlations. 
S 12 (q) is given by: 

5i 2 (g) = - > v - — ~ 2 ? s , 2 , r , (A9) 

,£41-^2(9)^(^7^) 

which is 0($i) and hence again generally acts as a perturbation. The explicit evaluation of 
the sums to calculate Sn(q) and Si 2 (q) is performed in an analogous manner to the single 
layer described by Bohm and Schweng£3 

We can note in passing that the success of the relatively simple H-RPA approach when 
applied to bilayer systems can be traced directly to Eqs. [A^ and |Aj|: the H-RPA includes 

31 



quite effectively the 0($i) contribution to the local field in a bilayer system. The STLS 
formalism, however, includes the inter-layer exchange-correlation and the perturbation of 
the intra-layer correlations by the presence of the second layer. 

We can rewrite Eq. [L^ for the local field factors Gij(q) for the case of equal densities in 
terms of dimensionless variables (Q = q/kp, K = k/k F ) as: 

g«(Q) = I / dK |Q||K| ^15 fe " ^ (IQ " K|)1 ' (A10) 

where S^(|Q — K|) = SV,-(|q — k|) and F^(|Q|) = i^-(|q|). Following Jonson in using the 
substitution C = Q — K, and with a = C/Q, we have as our final form:0 



Gv{Q) = z- !°° dC C [5^ - S'^C)] 
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^■(|Q|y / l + a 2 + 2acos0) acos0 
X / ^(IQI) y/l + a 2 + 2acoscl> ' 



Eqs. |A2| , [A3| and |A11| were solved self consistently using a method due originally to Ng 
to accelerate (ensure) convergence B For low number densities (i.e. r s > 2) it was also found 
that the initial value of Gu(q) has to be chosen carefully in order to ensure convergence; 



parameterised forms given in Ref. were used as the initial seeds. 
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FIGURES 




FIG. 1. The intra-layer and inter-layer Coulomb form factors Fn and F12, determined with 
quantum size effects taken into account (•). Fits using the functional forms given by Eq. ^ are in- 
dicated by the solid lines. The dashed lines give the form factors F^q) = 1 and F® 2 (q) = exp(— qd) 
for the strictly two-dimensional approximation. Note the logarithmic scale for F\2- 



36 




FIG. 2. Variation of (a) the intra-layer Gn(q) and (b) inter-layer Gu(q) local field factors 
with wavevector q, for 6 = 0.05 (solid line), 1.0 (dotted line) and 2.0 (dashed line); 9 = T e /Tp, 
T F = 78 K. 
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FIG. 3. Variation of the intra-layer and inter-layer local field factors G\\{q) and G\2{q) with 
electron temperature T e , for q = 0.2 hp. Note the scaling factor of 10 3 for Gu- 
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FIG. 4. Variation of the pair correlation functions gn(r) and g\2 (r), with in-plane separation 
r, for 6 = 0.05 (solid line), 1.0 (dotted line) and 2.0 (dashed line); 6 = T C /T F , T F = 78 K. 
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FIG. 5. A comparison of the T = local field corrections G\i(q) and Gu(q) determined in the 
ideal 2D approximation (dashed line) and with quantum size effects of the quasi-2DEGs taken into 
account (solid line). 
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FIG. 6. The pair correlation functions, gu(r) and gi2(r), as for Fig. ||. 
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FIG. 7. Theoretical Raman spectra for charge-density fluctuations (solid lines) and SPE 
(dashed lines) for T e = 10 K (0.13 Tp) and T e = 50 K (0.64 Tp). The damping parameter 
7 = 0.09 meV and the wavevector q = 1.6 x 10 5 cm -1 . The occupation factor (n(u>) + 1) in 



Eq. 11 has been omitted. Peaks due to the acoustic (AP) and optic (OP) plasmons are present in 
the charge-density fluctuation spectra. With increasing T e , the SPE continuum spreads to higher 
energies, overlapping the AP energy and leading to an asymmetric broadening of the AP peak, 
resulting from Landau damping. 
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FIG. 8. (a) Bulk PL signal for T e = 26 K and magnified (x4) view of the quantum well PL. 
The solid line indicates the fit described in the text. The inset shows the variation of the 2DEG 
temperature T e with the lattice temperature Tl. (b) The log of the quantum well PL intensity (•) 
for various temperatures T e (spectra are displaced vertically for clarity). Fits to the data, using 



Eq. 18 are indicated by the solid lines. 
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FIG. 9. (a) A spectrum measured for an electron temperature T e = 93 K is shown in the 
upper panel and the spectrum corrected for the Raman occupation factor {n(uj) + 1) is shown in 
the lower panel; q = 1.6 x 10 5 cm -1 . The points show the experimental spectrum; the solid lines 
indicate the fit to the spectrum (mainly obscured by the points due to the good fit) and the fit 
with the exponential backgrounds due to the bulk and hot PL removed. In the upper panel the 
bulk and hot PL contributions are indicated by the dashed lines, and the total PL contribution 
by the dotted line. The HWHMs (half- width at half maximum), y\ and 1/2, used to parameterise 
the peak widths, are indicated, (b) Corrected Raman spectra (with the background PL spectra 
subtracted and the (n(uj) + 1) occupation factor removed) of the acoustic (AP) and optic (AP) 
plasmons for various electron temperatures T c ; q = 1.6 X 10 5 cm" 1 . Note the broadening and 
developing asymmetry of the AP with increasing T c , resulting from Landau damping. 
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FIG. 10. Experimental (•) and theoretical (curves) results for the total width of the AP as 
a function of T c . Calculations have been performed within the RPA (dashed line), Hubbard ap- 
proximation (dotted line) and finite temperature STLS (solid line), (a) q = 1.16 x 10 5 cm -1 , (b) 
q = 1.6 x 10 5 cm -1 . 
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FIG. 12. The degree of asymmetry of the AP Raman peak, a = (y\ — j/2) / (ui + 2/2), for experi- 
ment (points) and theory, within the finite temperature STLS formalism (curves), for wavevectors 
q = 1.16 (o, dashed line) and q = 1.6 x 10 5 cm -1 (•, solid line). 
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FIG. 14. Experimental (•) and theroretical (curves) temperature (T e ) dependence for 
q = 1.6 X 10 5 cm -1 of (a) the OP energy wop; (b) the AP energy wap! (c) the AP damping 
Jwap = (yi + 2/2)- The theoretical curves have been determined for electron densities of (solid 
lines) 1.95, (dashed lines) 1.7 and (dotted lines) 1.5 x 10 11 cm -2 , using STLS local field factors 
given in Fig. ||. 
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FIG. 15. Experimental (•) and theroretical (curves) temperature (T e ) dependence for 
q = 1.6 x 10 5 cm -1 of (a) the OP energy loop', (b) the AP energy u ap; (c) the AP damping 
5to ap = (yi + U2)- The theoretical curves have been determined within the RPA using a Hubbard 
local field correction for a total electron density N\ + N2 = 3.9 x 10 11 cm~ 2 with N1/N2 = 1.0 
(solid lines), 1.9 (dashed lines), 2.3 (dotted lines) and 4.0 (dot-dashed lines). 
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FIG. 16. The temperature (T e ) dependence of the real part G' u (•) and imaginary part G'^ 
(o) of the intra-layer complex local field factor Gu(g,w), used as a fitting parameter to provide 
agreement between experiment and theory for the AP energy wap and the AP damping 5uj ap- 
The inter-layer local field factor is set to Gu = (a valid approximation given its relatively small 
magnitude), (a) q = 1.16 x 10 5 cm -1 , (b) q = 1.6 X 10 5 cm -1 . 
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